function PICP = calc_PICP(actual, lower, upper)
    within_interval = (actual >= lower) & (actual <= upper);
    PICP = mean(within_interval);
end